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We calculate the friction of fully mobile graphene flakes sliding on graphite. For incommensurately 
stacked flakes, we find a sudden and reversible increase in friction with load, in agreement with 
experimental observations. The transition from smooth sliding to stick-slip and the corresponding 
increase in friction is neither due to rotations to commensurate contact nor to dislocations but to 
a pinning caused by vertical distortions of edge atoms also when they are saturated by Hydrogen. 

This behavior should apply to all layered materials with strong in-plane bonding. 

PACS numbers: 68.35.Af,62.20.Qp,68.37.Ps,81.05.uf 


I. INTRODUCTION 

Since the discovery of graphenoi, there is a growing 
interest in materials similar to graphite, with strongly 
bonded 2D sheets and weak interplanar coupling^. This 
type of bonding makes these lamellar materials, and 
graphite in particular, also interesting as solid lubricants. 
Since the first Friction Force Microscope (FFM) measure¬ 
ments by Mate et al4 graphite has been a prototype sys¬ 
tem for nanotribology, the new research field that aims at 
understanding the fundamental origin of friction at the 
atomic scale. To model the friction in graphitic materi¬ 
als, it is important to go beyond the Tomlinson model 
of a point-like tip 4- ® and consider the tip as an extended 
contact. In fact, to account for the order of magnitude 
of the experimentally measured friction on graphite, it 
was proposed that a graphene flake was dragged along 
with the tip so that the friction between a graphene flake 
and graphite was measured instead of that between tip 
and graphite^ - — . For such an extended contact, the fric¬ 
tion strongly depends on the orientation of the flake with 
respect to the substrate. When the flake has the same 
orientation as the substrate, the contact is commensu¬ 
rate with high energy barriers to sliding and thus high 
friction. This is the most energetically favorable configu¬ 
ration, except near the edges of the substrate^. For an¬ 
gles in between commensurate situations, the flake is ap¬ 
proximately incommensurate and the potential barriers 
to sliding are averaged out. In this situation, an almost 
frictionless sliding was observed^ ’ 11 1 12 , a phenomenon of¬ 
ten referred to as super lubricity^ - —. 

Recently, extremely high speed superlubricity has been 
observed for micron sized graphene flakesii. However, 
for finite flakes, superlubricity is not necessarily a stable 
state. Further experiments and numerical simulations 


showed that graphene flakes of the order of 100 carbon 
atoms very often rotate to the commensurate orienta¬ 
tion with a large and irreversible increase of frictioi*i&. 
This finding was further supported by theoretical work, 
although stable orbits are predicted to exists. 

A sudden increase of friction was also measured for in¬ 
commensurate graphene flakes with increasing load- —, as 
shown in Fig. [I] Since the observed increase in friction 
was found to be reversible, it could not be explained by 
rotations to the commensurate state. Reversibility also 
rules out plastic deformations of the substrate. Models 
with rigid flakes cannot explain this increase since, for in¬ 
commensurate contacts, the corrugation remains too flat 
for the occurrence of stick-slip instabilities^^. Bonelli 
et alSi did simulations on non-rigid graphene flakes, but 
they only discussed the influence of load strongly limit¬ 
ing the deformations of the flake by means of very stiff 
springs (A'=2.5 eV/A 2 ). 

Under load, a breakdown of superlubricity can also oc¬ 
cur if strong in-plane distortions lead to local commensu- 
rability. This effect was reported by Kim and Falk^ for a 
model system of atoms with Lennard-Jones and harmonic 
interactions. They showed that the tip would adjust to 
the substrate with higher load or weaker harmonic inter¬ 
action between tip atoms. This adjustment led to local 
commensurability and consequently to the breakdown of 
superlubricity. Here we show that the reversible increase 
of friction with load shown in Fig. |T]is instead due to pin¬ 
ning of the edge atoms involving mostly vertical motion 
and very little in-plane strain in the flake. This mecha¬ 
nism seems to be specific of lamellar materials where the 
creation of defects or dislocations is energetically very 
unfavorable, contrary to the case of metal and rare gas 
islands on surfaces where the occurrence of dislocations 
dominates the diffusion^. 
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FIG. 1. FFM measurements of the load dependence of the 
friction force between a graphene flake on an HOPG graphite 
surface over a region of the graphite to which the flake is in¬ 
commensurate. Every data point was the average over 5 sep¬ 
arate measurements. The sudden increase of friction at loads 
in excess of ~ 40 nN was reproducible and fully reversible. 


II. MODEL 

We construct a model of a FFM experiment where 
a graphene flake made of N atoms is attached through 
springs to a tip that is moved on a graphite substrate. 
We allow the atoms of the flake to move in all directions 
whereas we keep the substrate atoms at their equilibrium 
positions in a flat hexagonal lattice at z = 0. The tip is 
modeled by attaching, in the x and y directions, each 
atom of the flake to springs at the positions of a rigid 
support flake of the same shape as shown in Fig. [5Ji. 
Alternatively, one could attach a spring to the center of 
mass, but this would not limit the rotations of the flake 
towards a commensurate contact, which is not the phys¬ 
ical situation we want to consider. We model the load as 
a constant force on each atom and report either the load 
per atom or the total load on the flake given as the sum 
of the load on all atoms. The interaction between the 
atoms in the flake is given by the REBO potentia l 24 ^ 5 
as implemented in the molecular dynamics (MD) code 
LAMMPS2&. The equilibrium interatomic distance d is 
1.3978 A which is close to the experimental value 1.42 A, 
giving a periodicity in the x direction a = a/3 d = 2.42 A. 
We describe the interlayer interactions in graphitic sys¬ 
tems by means of the the Kolmogorov-Crespi (KC) po¬ 
tential^. This combination has been shown to accurately 
reproduce the potential energy surface due to the sub¬ 
strate^ which is underestimated by the Lennard-Jones 
potential in AIREBO 24 . With this potential the interpla- 
nar distance in graphite is 3.34 A, with an energy gain of 
48 meV/atom and a difference of 15 meV/atom between 
A A and AB stacking. 

The forces acting on the flake atoms are given by 

Fife + ^ rebo + F spring H - F load (i) 

where the spring force on atom i is 

F spring,! = A (r ? ; I*/ (2) 


where Fj and F° are the in plane coordinates of the 
flake atom and its support point respectively, as indi¬ 
cated in Fig. [2ji. The spring constant K is taken to be 

16 meV/A 2 . 

The component in the pulling direction x of the total 
spring force, 


N 


f x = J2 f s 


i =1 


( 3 ) 


is often called the lateral force and its average over time 
gives the friction force Ff r i C . In fact, it can be shown that 
the average of the lateral force over a period of length a 
is 

AW 

^ F x >= Ff r i c = (4) 

a 

where AW is equal to the energy dissipated over a pe¬ 
riod of length o22. In our simulations, we use the first 
period as a transient and evaluate F/hc by averaging over 
the subsequent three periods. Notice that, in absence of 
interactions with the substrate, F x = 0 although each 
spring force F spr i ng f can be different from zero due to 
relaxation of the flake induced by the edge termination 
with respect to the fixed support points. 

The load on atom i is a constant force 


F loadti = -L/Nz (5) 

where L is the total load on the flake. 

Since the contact area of the flake has been estimated 
to be of the order of a hundred atomsSi^U, we choose flakes 
with hexagonal symmetry consisting of N = 54, 96 and 
150 atoms. We will mostly consider the case where the 
flake is rotated by 30 degrees with respect to the sub¬ 
strate to ensure incommensurability as shown in Fig.[2jr. 
This angle corresponds to an incommensurate contact 
for infinite lattices. By shifting the starting position of 
the incommensurate flake orthogonally to the pulling di¬ 
rection we examine the different scanlines indicated in 
Fig. CD). 

The speed of a FFM is so low ( ~ 1 — 1000 nrn/s) that, 
to a good approximation, the movement can be consid¬ 
ered static. We therefore follow the approach by Bonelli 
et ab 21 ' and use a quasistatic approach in which the sup¬ 
port is moved in steps along a given scanline and the flake 
is relaxed after each step by minimizing the force (Eq. 1) 
on each atom. We use the FIRE schemed, a damped 
dynamics algorithm, as implemented in LAMMPS2£. We 
use 200 minimizations per period, which gives a step size 
for the movement of the support of 0.012 A. The min¬ 
imization stops when the norm of the global force vec¬ 
tor becomes less than 3.12-10 -4 eV/A. We repeat this 
procedure for 4 periods to confirm periodicity and eval¬ 
uate Ff ric as described above. We found that the FIRE 
scheme is superior to conjugate gradients methods in sat¬ 
isfying the expected periodicity of the motion. 

In comparison with tight binding^, our model is more 
suitable to study the effect of load due to the longer cut¬ 
off. The 3.34 A interlayer distance in graphite is much 
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FIG. 2. (Color online) a) Model of the rigid tip and mobile 
flake. Each flake atom with coordinate r; is attached with a 
spring to a support point rj 1 . b) Flake with <f> = 30° on the 
substrate. The scanlines over which the center of mass of the 
support moves, from top to bottom: A, B, C, D. The period 
a is indicated. 
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FIG. 3. (Color online) Friction and friction coefficient [i as a 
function of load for commensurately oriented hexagonal flakes 
of different sizes as indicated. 


longer than the cut-off range used in the tight binding 
model (2.6 A). 

The quasistatic approach does not include dynamic 
and temperature effects. In some cases, we therefore 
compare results obtained by this method to MD simu¬ 
lations at constant temperature using a Langevin ther¬ 
mostat. We choose a damping parameter 7 -1 of 0.6 ps, 
a timestep of 1 fs, temperatures of 10 K and 300 K and 
move the support with a speed of v = 4.84 m/s. 


III. RESULTS 

In general, incommensurate contacts lead to much 
lower friction than commensurate ones. In Fig. [3] we 
show the friction as a function of load for Hakes with 
orientation <f> = 0, commensurate to the substrate. The 
friction linearly increases with load, resulting in a nearly 
constant friction coefficient fj, = Ff r i C /L = 0.03, which 
agrees well with the experimental value found for mi¬ 
croscale graphene^,. It is worth noting that /z slightly 
decreases with size as a result of the fact that atoms at 
the edges can reach deeper minima and contribute more 
to the friction. 

In experiments, the point of zero friction would occur 
at negative values of load due to the attractive van der 
Waals forces between tip and substrate. As in our sim¬ 
ulations only the contact area of the tip and substrate 
is modeled, the van der Waals forces are very small and 
already at L = 0 the friction is nearly zero. 

The effect of the edges becomes more dramatic for 
incommensurate contacts. In this case, the friction is 
drastically different from the one just discussed for com¬ 
mensurate cases. In Fig. [4] we show the variation of 
the friction force with increasing load for the 54, 96 and 
150 atom flakes and the A,B,C,D scanlines. Smaller, 24 
atom, flakes are not included because they rotate to the 
commensurate orientation already at low load, due to the 
smaller moment of inertia. We see that at low load, the 
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FIG. 4. (Color online) Friction as a function of load for several 
scanlines and flake sizes of a flake incommensurately oriented 
with 4> = 30° (left). The line labeled ’H’ shows the friction of 
a C 54 H 16 flake along scanline D. 


friction is indeed almost vanishing as expected for truly 
incommensurate case s 16 : 33 . 

At higher loads, not only is the friction at least one 
order of magnitude lower than for commensurate flakes 
but the dependence on load is also much more complex. 
In most cases, the friction suddenly increases from a cer¬ 
tain threshold load onwards in a way similar to the ex¬ 
perimental situation shown in Fig. [lji. The increase in 
friction is stronger, and starts at a slightly lower load, 
for smaller flakes. For the 96 atom flakes, the increase 
starts at L/N ~ 0.6 nN or L k 60 nN, which is close 
to the experimental valued. The large increase in fric¬ 
tion with load corresponds to the transition from smooth 
motion to stick-slip behavior shown in Fig. [5] For an in¬ 
commensurate contact, stick-slip motion is very unusual 
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FIG. 5. (Color online) Lateral force as a function of the sup¬ 
port position for a low (blue, gray) and a high (red,black) 
load for the A,B,C,D scanlines (top to bottom) and different 
sizes (left to right). The total load L is indicated in the figure. 
Notice the change from smooth to stick-slip behavior. 


and would only be expected for very small flakes which, 
strictly speaking, are not incommensurate with the sub¬ 
strate^. 

Fig. S]also shows that the behavior with load is not the 
same for all scanlines and flake sizes. First, we examine 
in detail the simplest case to explain the behavior with 
load and the effect of deformations and thereafter we will 
discuss the general features for all scanlines and flake 
sizes. 

The simplest case is presented by the 54 atom flake 
along scanline D because the flake remains at 30 degrees 
and its center of mass does not deviate from the scan¬ 
line, as shown in Fig. [6] We will focus on the role of 
deformations of the flake which are usually neglected. 

A priori, it is not obvious whether to expect higher or 
lower friction by considering a non-rigid flake instead of 
a rigid one. While barriers are easier to circumvent by 
deforming in addition to shifting and rotating, the same 
freedom could be used to find deeper minima. Therefore 
we examine the energy profile before looking at the fric¬ 
tion. In Fig. [7}i we show the potential energy due to 
the substrate calculated for a rigid flake kept at a fixed 
height, corresponding to the average height at the given 
load. We see that there is a smooth minimum at a/2, 
that becomes deeper with load. To show the importance 
of the deformation of the flake, we compare this result to 
the energy profile for a deformable flake. To obtain this, 
instead of pulling the flake through the coupling to the 







x(A) 


tip 


(A) 


FIG. 6. (Color online) Two periods of the trajectory of the 
center of mass (left) and rotation angle versus the tip position 
(right). Simulations are for the 54 atom flake with L = 80 nN 
or L/IV=1.48 nN (top), for an 96 atom flake with L=140 nN or 
L/N=1A6 nN (middle) and for a 150 atom flake with L=220 
nN or L/N= 1.47 nN (bottom). Since the flakes are coupled 
to the support with springs, the center of mass can deviate 
from the center of mass of the support (shown in black). 


support, we shift the support with the attached flake and 
minimize the energy always from the same starting flake 
configuration. In Figs. [7jr and [7J:: it can be seen that the 
energy profile E tot obtained in this way is very different 
from the one for a rigid flake (Fig. [TJa,). This is mainly 
due to the contribution to the energy associated with 
load, El = L ■ zcm- As load increases from 20 nN (Fig. 
(TJd) to 40 nN (Fig. [7J:) the minima in the energy profile 
of the relaxed flake become deeper and, what is more im¬ 
portant, a sharp barrier at Xu p = a appears in between 
at high load. This barrier gives rise to a discontinuity in 
the difference between the center of mass of the flake and 
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the support shown in Fig. [7J1. When we pull the flake 
using the quasistatic method described in section II, the 
sharp barrier at high loads makes the flake stick instead 
of smoothly following the support. The change of behav¬ 
ior from continuous to stick-slip is evident in the total 
energy as a function of tip position in Fig. [7J:. At a load 
of 20 nN the barriers are not curved steeply enough to 
pin the flake, resulting in smooth movement, whereas at 
a load of 40 nN the flake remains pinned in the minima 
shown in Fig. [TJc and the motion becomes discontinu¬ 
ous. As load increases further, a second slip emerges as 
also the second transition between the two minima is no 
longer barrier less. The dominant contribution to these 
energy profiles is the load, as a consequence of atoms 
being closer or farther away from the substrate. 

When we compare the average distance to the surface 
of the edge atoms, defined as the atoms with two neigh¬ 
bors, to that of the inner atoms (Figs. [Tf and [TJ*), we 
see that the edge atoms are much closer to the surface 
and therefore contribute most to the load energy. While 
the inner atoms move nearly continuously also at higher 
loads, the edge atoms move discontinuously. The two 
configurations corresponding to the two minima in the 
period a at L = 40 nN are shown in Fig. [8] One can 
see that edge atoms on either side are much closer to 
the surface and the discontinuity of the motion corre¬ 
sponds in this case to a tilting of the flake in going from 
Xtip = 0.315a to 0.685a and vice versa. 

All examined flake sizes and scanlines which show a 
transition from smooth to stick-slip motion as a function 
of load, show a corresponding jump in the z coordinates 
of the edge atoms. As the edge to surface ratio decreases 
with flake size, edge effects become less important for 
larger flakes. We indeed observe that the friction in¬ 
creases more, and from a lower load per atom, for smaller 
flakes, as shown in Fig. [U That the increase in friction 
at incommensurate orientations is mostly due to edge ef¬ 
fects makes friction very dependent on the details of the 
energy landscape and makes it difficult to draw general 
trends as a function of scanline, number of particles, and 
orientation with respect to the substrate. Depending on 
the scanline, there can be multiple minima or only one 
per period and the barriers between them can be high or 
low. For instance for the 96 atom flake, a marked increase 
of friction is observed for scanlines C and D whereas no 
or only a very small increase is found for scanlines A and 
B respectively. The increase also depends on the orienta¬ 
tion of the flake. For less symmetric situations, like flakes 
with (f> = 15° and 4> = 25°, we do observe an increase in 
friction with load for all scanlines. 

For the 96 atom flake along scanline C, we have also 
performed MD simulations. In Fig. [9] we see that at 
10 K the flake displays the same stick-slip behavior found 
in the quasistatic approach. At room temperature how¬ 
ever, the stick-slip motion is masked by large fluctuations 
of the lateral force, resulting in negligible friction. As 
the speed of nearly 5 m/s is several orders of magnitude 
higher than in experiments, these simulations are likely 
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FIG. 7. (Color online) 54 atom flake along scanline D: (a) En¬ 
ergy as a function of the tip position for a rigid flake at a fixed 
height z = 2.97 A and 2.80 A for L = 20 nN and L = 40 nN 
respectively, (b) Different contributions to the energy as a 
function of the tip position for a deformable flake which is 
relaxed from its ideal configuration for each tip position, for 
L = 20 nN. (c) Same as in (b) for L =40 nN. (d) Distance 
between the center of mass of the relaxed flake and the center 
of mass of the tip. (e) Variation of the total energy from that 
at Xtip = 0 for a pulled flake for several loads, (f) The average 
distance to the substrate for the edge and non-edge atoms for 
L = 20 nN. (g) Same as in (f) for L = 40 nN. 
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FIG. 8. (Color online) Two minimum energy configurations 
in the trajectory of the 54 atom flake along scanline D for 
L = 40 nN. The color shows A z = Zi — zcm where cm refers 
to the center of mass of the flake: red is closer to the substrate, 
blue is farther away. The outer atoms are clearly more mobile. 


the flake and substrate atoms through the REBO poten¬ 
tial and neglect interaction with the tip. The effect of 
Hydrogen on friction as a function of load for this flake 
along scanline D is shown by the line labeled ’H’ in fig¬ 
ure 4. The behavior of friction is qualitatively the same 
but the increase starts at a higher load when Hydrogen 
is present. This fact might have been expected because 
the in-plane bonds are better preserved up to the edges, 
making the edges of the flake less flexible. 

In addition to the effects studied here, elastic deforma¬ 
tions of the substrate could also add to the total friction. 
The latter effect decreases with the number of layers of 
the substrat o 35 ’ 36 . As we model a bulk graphite sub¬ 
strate, we have not included this in our model. However, 
elastic deformations could significantly contribute to the 
friction on few-layer graphene substrates. 


to underestimate the friction. The motion of the flake 
in the stick-slip regime is related to a very symmetric 
configuration of the flake with the six corner atoms lock¬ 
ing into favorable positions of the substrate potential as 
shown in the right panel of Fig. [9] 



position for a 96 atom flake along scanline C for L = 100 nN 
or L/N ~ 1 nN calculated by MD at 10 K. The constant 
lines give the resulting friction at 10 K and the one at 300 K 
( F x not shown). The asterisk corresponds to the minimum 
energy configuration shown in the right figure. The red (dark) 
color indicates atoms closer to the substrate as in Fig. [8]with 
minimum at Az = —0.165 A. 


IV. SUMMARY AND CONCLUSION 

In summary, we have studied the motion and friction 
of mobile and flexible graphene flakes moving at incom¬ 
mensurate orientations on a graphite substrate. In agree¬ 
ment with FFM experimental results, we have found that 
the superlubric behavior at low loads evolves to a fric¬ 
tional stick-slip motion at high loads. This change is 
reversible since it is not due to rotations to commensu¬ 
rate orientations but rather to a locking of the flakes as 
a result of vertical motion of the edge atoms, also when 
H-saturation is considered. Interestingly no dislocations 
appear in the flake also under high load, contrary to the 
typical behavior of diffusion for metals and rare gas is¬ 
lands on surfaces. The strong in-plane bonding of layered 
materials like graphite makes that the crystalline struc¬ 
ture is preserved while vertical distortions at the edges 
are energetically favorable. This feature might explain 
the good and persisting lubricant properties of layered 
materials. 
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